Equilibrium and nonequilibrium fluctuations at the interface between two fluid phases. 
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We have performed small-angle light-scattering measurements of the static structure factor of a 
critical binary mixture undergoing diffusive partial remixing. An uncommon scattering geometry 
integrates the structure factor over the sample thickness, allowing different regions of the concentra- 
tion profile to be probed simultaneously. Our experiment shows the existence of interface capillary 
waves throughout the macroscopic evolution to an equilibrium interface, and allows to derive the 
time evolution of surface tension. Interfacial properties are shown to attain their equilibrium values 
quickly compared to the system's macroscopic equilibration time. 
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I. INTRODUCTION 

The fluctuations at the interface between two fluid phases at thermodynamic equilibrium have been studied very 
extensively starting from the beginning of this century |l| , and particular investigation has concerned those at the 
interfaces of critical fluids [^-^ . Although the features of equilibrium interfacial fluctuations are now relatively well 
known, the behavior of an interface under nonequilibrium conditions is still not well understood. Many experiments 
have been performed to detect an effective nonequilibrium surface tension in miscible fluids [^-§1 but the results 
obtained are mainly qualitative. 

In this paper we will present experimental results about the behavior of interfacial fluctuations during the diffusive 
remixing of partially miscible phases. In this system the interface between two fluid phases is being crossed by a 
macroscopic mass flow. It is well known that the fluctuations in the equilibrium states before and after the diffusive 
partial remixing are controlled by surface tension and gravity. It is not clear, however, what happens to the interface 
during the remixing: does the interface temporarily dissolve? Is there a surface tension during the transient, and 
how is it related to the evolving macroscopic state? We try to address these problems by means of low-angle light 
scattering measurements of the correlation function of fluctuations. 

The sample considered is a near critical binary mixture kept below its critical temperature Tc, so that it is macro- 
scopically separated into two bulk phases by a sharp horizontal interface. The diffusive remixing is started by raising 
^ , the temperature to a value closer to Tc, but still below it. 

• We report data at equilibrium showing the expected interface capillary waves, and data obtained out of equilibrium 
Ch ', during partial diffusive remixing, showing for the first time that capillary waves at the interface are still present 
O ' and coexist with "giant" nonequilibrium fluctuations in the bulk. By combining time- resolved measurements with 
predictions for the light scattered by the fluctuations, we are in the position to derive data for the time evolution of 
the nonequilibrium surface tension. These data show for the first time that during the diffusive remixing the surface 
tension attains almost instantaneously its final equilibrium value, and this is consistent with a fast rearrangement of 
the concentration profile in the neighborhood of the interface. 
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1 



II. EXPERIMENT 



Traditionally, interfacial fluctuations have been studied by means of dynamic surface light scattering techniques [ p^ . 
These techniques allow to determine the power spectrum of the light scattered from the excitations, and have been 
used very extensively to characterize the equilibrium properties of interfacial fluctuations in simple fluids and binary 
mixtures . Surface light scattering is usually performed by sending a probe beam on the interface at the angle 

of total reflection, scattered light being collected around the specular reflection angle and partially recombined with 
the main beam, so that a heterodyne signal is obtained. However, dynamic scattering techniques are not very well 
suited to study time-dependent nonequilibrium processes, as the time needed to accumulate an adequate statistics of 
the fluctuations is often much longer than the time related to changes in the macroscopic state of the fluid. To bypass 
this problem we have used a unique low-angle static light scattering setup. Although the fluctuations' timescales 
are not accessible with this instrument, it allows to determine in a fraction of a second the static light scattered at 
31 wavevectors distributed within a two decades range, making it a very useful tool to study fluctuations around 
a time-dependent macroscopic state. This instrument, described in detail elsewhere | p^ , p^ , typically investigates a 
wavevector range corresponding to lOOcm"^ < 9 < lOOOOcm"^. 

Our setup is configured to detect the transmitted static scattered intensity, the probe beam being sent vertically at 
normal incidence. In this way light is scattered both at the interface and in the bulk layers above and below it. It can 
be easily shown that the structure factor in the transmission geometry is proportional to the usual one in reflection. 

A critical binary mixture is an ideal sample to study the nonequilibrium fluctuations during diffusion at an interface. 
Although other partially miscible fluids could be used to perform this experiment, the use of a critical mixture allows 
to tune the timescale of the macroscopic diffusion process by adjusting the temperature difference from the critical 
point. Moreover experimental runs can be iterated simply by cycling the temperature. 

The sample is a 4.5 mm thick horizontal layer of the binary mixture aniline-cyclohexane prepared at its critical 
consolution concentration (c=0.47 w/w aniline). Its critical temperature Tc is about 30°C, and was determined to 
within ±0.01°C before each experimental run. While slowly decreasing the temperature of the single phase above T^,, 
Tc was taken as the temperature where a sudden increase of turbidity was observed. 

The light scattering cell is a modification of the Rayleigh-Benard one already used to investigate fluctuations in a 



thermal diffusion process 
between two massive sapp 



configured to keep the sample at a uniform temperature. The mixture is sandwiched 
)hire windows whose temperature control, achieved with Peltier plates, is as good as 3mK 
over a period of one week, temperature differences between the two plates being kept to about 2mK. 
A typical measurement sequence involves the following procedure. An optical background is recorded with the mixture 
in its one phase region at 5.5 K above Tc, where the bulk fluctuations amplitude is many orders of magnitude smaller 
than the one of equilibrium and nonequilibrium fluctuations in our experiment. This optical background mostly 
contains contributions due to dust and imperfections of the optical elements, and it is subtracted to all subsequent 
measurements. The system is then let phase separate at 3.5K below Tc, where the concentration difference between 
the two bulk phases is Ac ~ 0.5 |0. Great care is dedicated to eliminate wetting drops at the optical windows. After 
a few hours the intensity distribution scattered by the system at thermodynamic equilibrium is recorded, scattered 
light being mostly due to the capillary waves at the interface. The nonequilibrium process is then started by suddenly 
increasing the temperature at 0.1°C below Tc. The concentration difference between the phases has to readjust by 
means of a diffusive process. The two macroscopic phases at this temperature are not completely miscible. At equi- 
librium they will have concentration difference of Ac « 0.15 [|lj] and they will be separated by a new interface. The 
scattered intensity distribution is recorded during this transient, until the system reaches thermodynamic equilibrium 
after about 24 hours. It is well known that the interfacial fluctuations at the equilibrium states preceding and following 
the diffusion process are overdamped capillary waves JTof , characterized by a power spectrum which exhibits a 
gravitational stabilization at small wavevectors. The main problem we want to address is how the fluctuations behave 
during the transient between the equilibrium states. 
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III. DISCUSSION 



The evolution of the scattered mtensity distributions can be roughly divided into two stages. In the first stage, 
represented in Fig. 1, the intensity distribution, initially due to the equilibrium interface excitations, increases and 
develops a bump, which represents the appearance of a typical lengthscale. The initial data set, being the least 
intense, is particularly affected by the subtraction of the optical background. The time elapsed during this stage 
roughly corresponds to the thermal time needed to increase the temperature of the sample (about 200 s). In the 
second stage, represented in Fig. 2, the intensity distribution decreases, until eventually the bump disappears, and 
the q~'^ equilibrium power spectrum of capillary waves is recovered. Notice that light scattered at large wavevectors 
does not change with time during this stage. We will see that this is related to the readjustment of the concentration 
profile across the interface having taken place. 

To analyze our data we assume that the sample may be depicted as the superposition of two thick bulk layers 
separated by a thin interface layer, and that these layers scatter light independently from each other. Therefore we 
are considering two sources of scattering, the interface fiuctuations scattering Imtil) and the fluctuations in the bulk 
phases scattering Ihuikio), so that the total intensity distribution is: 



I{q) = Iint{q) + huik{q), 



(1) 



where, as outlined in the appendix, 

dljnt (q, t) 



dn 



(27r)2 yOcJ p(3g 



-kuT- 



qcap 



(2) 



and 



dlbuik{tl,t) 
dQ. 



n^K^ (dn 



(27r)2 \dcj pl3g 



1 , ^ Acbuifc 
-kbT ■ 



(3) 



In Eqs. (|^) and (^, is the wavevector of light in vacuum, n, c, p are the mixture's refraction index, weight-fraction 
concentration and density, respectively. /? is and g is the gravity acceleration. 

Aci„4 is the concentration difference across the interface and Ac;,tiifc is the total sample concentration difference minus 
Acint, that is the concentration difference that falls in the bulk phases. The roUoff wavevectors qcap and qro, given by 



qcap 



^P9 



(4) 



and 



2vD 



(5) 



characterize the onset of gravitational stabilization at large lengthscales of capillary and bulk fluctuations, respectively 
[0,^. In Eq. (^) dzC represents the largest concentration gradient in the bulk phases at a certain time ||l5| . 

We use Eqs. (0)-(||) to fit the experimental data shown in Figs. 1-2. In order to limit the number of fitting 
parameters we take advantage of the fact that the concentration difference Ac across the whole sample does not 
change during the so called free-diffusive regime [|l^,^ , as the diffusive remixing initially involves only layers of fluid 
close to the interface. Therefore Aci„t and Acbuifc are related by 



Aci„t(t = 0) = /S.Cint{t) + Acbulk{t) 



t < r„ 



(6) 



^jJ-Q is the time required for diffusion to occur over the sample height s. Tmacro is of the order 



where Tmacro 

of 7000s for our sample, by assuming D=6 • 10~^^^ (this is the equilibrium value at 3K below Tc, see and 
references therein). By imposing the reference value Aci„i(t = 0) = 0.5 |lj] and by fitting the experimental data 
using Eqs. (0)-(^, we are able to determine three parameters: the concentration difference across the interface and 
the roUoff wavevectors qcap and qro ■ 
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FIG. 1. Normalized scattered intensity plotted vs. scattered wavevector q at different times. The initial dataset (v) is the 
intensity scattered by the initial equilibrium interface. Time is measured from the start of the partial remixing process. The 
figure shows the early stages of the remixing process, when the scattered intensity is increasing with time. The solid lines 
represent the best fit of the experimental data with Eqs. (h 




FIG. 2. Normalized scattered intensity plotted vs. scattered wavevector q after the onset the remixing process. The scattered 
intensity is decreasing with time. The least intense scattering (y) is due to the final equilibrium interface. The solid lines 
represent the best fit of the experimental data with Eqs. (h 
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Results for Aci„t are presented in Fig. 3(a) as a function of time. 




FIG. 3. Time evolution of the concentration difference across the interface, (a), and of the interfacial surface tension, (b). 
Circles represent the experimental results obtained by fitting the scattered intensity distributions with Eqs. (|l]-|). Crosses 
mark the reference equilibrium values in the initial and final states. The triangle in panel (b) corresponds to the surface 
tension measured in the initial equilibrium state. Vertical bars on the experimental results represent the estimated fitting and 
sistematic errors. 



Initially Aci„t has its equilibrium value Aci„4 = 0.5. After about 40 seconds from the temperature increase it begins 
to drop, and decreases for about 300s, finally stabilizing to the constant value Acmj = 0.3. This is roughly a factor 
of two larger than the reference equilibrium value Q . This discrepancy is due to the difficulty (both theoretical and 
in the fitting procedure) of clearly ascribing the measured scattered intensity to either interface or bulk phases, and 
an estimate of our rather large fitting and systematic error is shown with the error bars in Fig. 3. This error is such 
that our data cannot be considered fully quantitative, but it does not affect the features we discuss. 
By combining the results for Aci„f and qcap, since from Eq. (^ 

P(3Acm9 

= 5 , (7) 

^cap 

we are in the unique position to obtain the time evolution of the interfacial surface tension during the nonequilibrium 
process. Experimental results for the surface tension are shown in Fig. 3(b), which represents the main accom- 
plishment of this work. The two crosses mark the value of the equilibrium surface tension at the initial and final 
temperature, extrapolated from the reference data from Atack and Rice [|l^. The first data point represents the 
equilibrium surface tension measured with our light scattering setup. The agreement with the reference value is good. 
After the diffusion process is started the surface tension drops about two orders of magnitude, until after about 300s 
it stabilizes to a constant value. The asymptotic value of the surface tension is about a factor of two larger than the 
reference value, a good result considering the wide range of values spanned. Although these results are only partially 
quantitative. Fig. 3 unambiguously shows for the first time that the properties of the nonequilibrium interface rapidly 
attain their equilibrium values. This equilibration time is very small compared with the one associated to readjust- 
ments of the bulk phases (which corresponds to about one day) , and it is comparable to the time needed to increase 
the temperature of the sample. Notice that the surface tension evolution docs not show the initial delay seen on the 
Acint evolution. This indicates that the surface tension is probably following the local temperature almost instantly, 
whereas Aci„t does not change until diffusion has occurred over the fluctuations' characteristic lenghtscales. With 
the diffusion coefficient given above this time is about 30s for the smallest wavevectors observed. 



We are now in the position to comment the fast growth, shown in Fig. 1, of the scattered intensity distributions at 
intermediate and large wavevectors. Soon after the diffusion process is started the concentration difference across the 
interface decreases at its equilibrium value. According to the concentration conservation Eq. M), a strong concentration 
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difference Acbuik is rapidly created in the bulk phases, and a large concentration gradient quickly grows near the 
interface. This gives rise to velocity-induced concentration fluctuations described by Eq. ^ [|15|,|6|. The growth of 
Acbuikit) during the free diffusive regime is shown by the circles in Fig. 4, and it mirrors the results for Aci„t in 
Fig. 3(a). When enough time has passed for diffusing particles to reach the macroscopic boundaries (about 5000s 
as from Fig. 4), the sample enters the restricted diffusion regime, where the total concentration difference across the 
sample begins to change, and Eq. (^ does not hold any more. However, according to Fig. 3, we can now assume 
that the interfacial parameters have attained their asymptotic values. In this way we can fit the scattered intensity 
distributions to determine the concentration difference across the bulk phases. Results for Acbuik obtained in this 
way are shown by the squares in Fig. 4. 



0.25 



<1 




FIG. 4. Time evolution of the concentration difference across the bulk phases. During the initial free-diffusion stage (circles) 
the total concentration across the sample is conserved (see Eq. ^) and the data mirror those presented in Fig. 3(a) for 
the concentration difference across the interface. Later on the presence of boundaries is felt, and the sample enters the 
restricted-diffusion stage (squares). In agreement with the results presented in Fig. 3, we have assumed that during this stage 
the interfacial properties have already attained their equilibrium values. 

The evolution of this process is very similar to that already observed during the free diffusion of completely miscible 
phases [ p"5p^ . From the fitting of the nonequilibrium bulk data we also determine the rolloff wavevector g^o, which 
corresponds to the bump observed in the light scattering data in Figures 1-2. As outlined above and thoroughly 
described in Ref . , fluctuations at wavevectors smaller than g^o are stabilized by gravity, which frustrates the 
divergence at small wavevectors. The rolloff wavevector is plotted in Fig. 5 as a function of time. 
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FIG. 5. Time evolution of the rolloff wavevector g^o- 

The results displayed in Fig. 5 show that the variation of q^o roughly corresponds to a factor of 1.4. g^o is mostly 
determined by the layers of fluid where the concentration gradient is largest ||l^ , that is the bulk layers close to the 
interface. From Eq. we can estimate that the variation of the concentration gradient close to the interface is 
rather small, roughly corresponding to a factor four. This behavior suggests that the bulk concentration gradient is 
pinned to the interface concentration profile. 

IV. CONCLUSIONS 

We have reported light-scattering measurements of the correlation function of fluctuations at the interface and in 
the bulk phases of a critical binary mixture undergoing diffusive partial remixing. From the intensity distributions it 
is clear that even during remixing a sharp interface separates the two bulk phases, and that this interface is roughened 
by capillary waves similar to those at equilibrium. Our data have been analyzed, yielding qualitative information 
on the concentration profile of the system. In particular our main results are the time evolution of the interface 
concentration difference and the surface tension, showing how these parameters relax to the final equilibrium values 
quickly compared to the system's macroscopic equilibration time. 
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APPENDIX: 

In this appendix we shall outline how one may derive the structure factor of fluctuations in an out of equilibrium 
binary mixture where a concentration gradient and a diffusive concentration flux are present. We consider a fluid 
described hy a. z dependent concentration gradient c{z). Layer by layer, surfaces of uniform concentration can be 
defined. From a macroscopic point of view these surfaces are horizontal planes. In the presence of fluctuations the 
surfaces are corrugated, due to the motion of parcels of fluid in the vertical direction. As customary, we will indicate 
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by hz{x,y,t) the vertical displacement of the surfaces from their mean position at z. 

In such a system light scattered with scattering vector q is proportional to the mean square amplitude of the roughness 
mode q, (|/iq(f)|^). We shall show how this quantity may be calculated for the two limiting cases of a sharp interface 
and a linear concentration gradient, extending the fluctuating hydrodynamics treatment applied in 
Throughout this paper we are dealing with small fluid velocities and overdamped motion, and equations will be 
approximated accordingly, see for example the discussion in [0. We will consider the scattering of light by long- 
wavelength fluctuations, which scatter light mostly at small scattering angles in the forward direction. 
We shall suppose fluctuations to be generated independently at different heights, and shall first consider a single 
fluctuation generated in a layer dz at z. Then, for every wavevector q, one has that hz (q) will be maximum for z = z 
where the fluctuation was generated, and exponentially decreasing with distance from z like hz = hz exp(— (j|z — z\). 
This is a consequence of the hydrodynamic equations for a viscous fluid in the case of overdamped motion pO| . 
The hydrodynamic equations describing motion in an incompressible viscous mixture can be linearized for small 
velocities H as: 



div V = 



(Al) 



d-\r{x,y,z,t) 
dt 



= i/V^v V P{x,y,z,t) 



(A2) 



dc_ 

'dt 



P 



(A3) 



where v is the fluctuating velocity, c the fluctuating concentration, v is the kinematic viscosity, P{x,y,z,t) is the 
pressure, (Vc)t i s th e concentration gradient averaged over typical fluctuations' timescales and j is the mass diffusion 
flux. Equation (Al) is the usual mass conservation equation for an incompressible fluid, Eq. (A2) is the linearized 
Navier-Stokes equation, and Eq. (A3) is the convection-diffusion equation. The macroscopic diffusive flux j is pro- 
portional to V/i, where fi is the chemical potential of the mixture; it is non-zero where fi is not constant. 
Equiconcentration surfaces are corrugated by thermal velocity fluctuations in the vertical direction, and these are 
described by the z component of Eq. ( A2 ) . 

We change variable from c to h in Eq. (A3), and simplify the equation by assuming that concentration diffusion obeys 
Fick's law ^ = DV^h, where D is the "mutual diffusion" coefficient. 

From Eq. (A2) we take the z component of the velocity, Vz — -^jh + DV^h, and substitute it into equation (|A2| ), 
obtaining, after Fourier transform in x,y and t, an equation of motion for non-equilibrium fluctuations: 



-iui [-iu! hq^^(z) + 2Dq hq^^{z)] = 
2i^q'^ [-iujhq,^{z) + 2D q'^ hq^^{z)] - ^qPq^^{z) + ^qSq,uj, 



(A4) 



where the continuity of tangential stresses on the fluctuating surface has been imposed and a stocastic force term Sq 
has been added to describe the onset of thermal spontaneous velocity fluctuations. The correlation function of this 



stocastic force is assumed to be (see Q) ( 



) 



pA 



as in equilibrium. 



The pressure P in Eq. (A4) may be exerted by the external gravity force and by internal capillary forces and is 
measured against the average fluid pressure. Because of Eqs. (Al)-(A2), in deriving Eq. ( A4) we have imposed that 
the pressure P{z) induced by the fluctuation in layer z decays exponentially with distance from the fluctuation layer 
to the average pressure, like Pq{z) — Pq{z) exp {— q\z — z\). This pressure term depends on the local concentration 
profile. It can be written explicitly for fluctuations involving a "sharp interface" profile and for those in the bulk 
phases, by considering the gravity and capillary forces on the fluctuations. We shall consider a fluctuation to be a 
sharp-interface fluctuation if z is close to the interface position Zi. Instead, if dzc{z) may be considered linear within 
a range of ±i around z, the fluctuation will be considered a bulk fluctuation. Our treatment is approximate in that 
we are supposing that every fluctuation falls within one of these cases. 
Interface fluctuations and bulk fluctuations give rise to different dynamics. 

In the case of a sharp interface, the pressure acting on a fluid element at height z — Zi, when a fluctuation occurs at 
Zi bending the interface, is: 



Pqi^i) = ^ [Apghq{z 



(A5) 
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Substituting this pressure in Eq. ( A4 ) one may, with algebra similar to that in , calculate the spectrum of interface 
fluctuations: 



2 a q + a- q + i p I' D q 

l^.l^ = (IMq)l')* , ""'^ (A6) 

^2 I / Apgq + aq^+4pi^Dq->- 
\ A pu q^ 



where the static term is: 



A/Xpg + crg^ + ^pvDq-^ 

At equilibrium, when the mutual diffusion coefficient = 0, the spectrum linewidth is the same as that of the 
well-known equilibrium interface H,|T^: gravity dominates at small wavevectors and capillary forces at large ones. 
However, during the nonequilibrium process, diffusion is effective in relaxing large-wavevector fluctuations. At these 
wavevectors the linewidth becomes the usual diffusive Dq^ one, and this diffusive contribution is also apparent in the 
static structure factor where a new term proportional to is present in the denominator. With our setup it was not 
possible to study the structure factor at large enough wavevectors to check the result and, as far as we know, this 
feature has never been observed. 

The other case is that of a bulk phase where the concentration gradient is approximately linear. Here the pressure 
opposing a fluid element at height z when a fluctuation occurs bending the layer at height z is: 

^q(^) = \ [d.pgh^iz) + dzaq^h^iz)] -, (A8) 
z q 

where dz stands for 

From mean field theories of binary systems |l^j2^ it follows that §f oc ) ■ If the concentration gradient is 
small, the capillary term becom es n egligible compared to the gravitational one, and we shall drop it from now on in 



1 kuT 



this case. As before, from Eq. (A4) one may calculate the correlation function of bulk fluctuations, recovering the 



spectrum calculated and commented in Ref. llC 



2 2dzpg + 4pyDq'^ 

\h^^{zt = (IMq, . ''"'^ m 

2d^pq + ApiyDq* 
4: p u q^ 



where the static term is 



Although these bulk fluctuations have the same origin as capillary waves, namely velocity fluctuations parallel to the 
concentration gradient, the transition from a sharp to a diffuse interface radically modifies both the static and the 
dynamic structure factor of the fluctuations. 

We shall now outline how one may evaluate the intensity of light scattered by these fluctuations. Suppose a plane 
wave having wavevector Kq in vacuum and intensity /q is propagating in the vertical direction in a sample having 
index of refraction n = n{z). The roughness of the equiconcentration surfaces due to a fluctuation at z introduces 
a phase dependancy on {x,y). One easily sees that the resulting scattered intensity per solid angle is — 
Iq ^(2^)" l^^i'i' where A$(q, t) is the optical path variation induced by the fluctuation given by 



A$(q,t) = ^ y dzhz{q,t) dzc{z). (All) 



This integral is similar to that required to calculate the pressures of Eqs. (A5) and (AS) and, depending on the 
system's local concentration profile at z, it can be easily approximated in the same two limiting cases considered 
above, sharp interface or concentration gradient in bulk phase. 

Scattering from the whole sample is given by integration over the sample thickness of the scattering due to fluctua- 
tions arising in a single layer, weighed with the probability of being in that layer and with the density of fluctuations. 
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so that one integrates the single fluctuation intensity in dz^. 

If the sample comprises both bulk regions and a sharp interface having a concentration difference Aci„i, the total 
scattered light is then: 



dl. 



q.t _ ^ n^K^ fdn 



d'^ 



(27r)2 \dc 
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I dz — 

hulk Q 


dzc{z) 




\hMM 




Z- 





(A12) 



We have used Eq. (A12) to fit our data after approximating the bulk phase integral by considering that most of the 
scattering is due to the layers with the greatest dzc{z). 
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